The g factor of an electron or muon bound by an arbitrary central potential 
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We consider the g factor of a spin- 1/2 particle (electron or muon) bound by an arbitrary central 
field. We present an approach which allows one to express the relativistic g factor in terms of the 
binding energy. We derive the general expression for the correction to the g factor caused by a 
deviation of the central potential from the Coulomb one. As the application of this method, we 
consider the corrections to the g factor due to the finite nuclear size, including vacuum polarization 
radiative correction. The effect of the anomalous magnetic moment is also taken into account. 

PACS numbers: 24.80.+y, 25.30.Bf, 21.10.Gv 



I. INTRODUCTION 

Study of energy levels of atomic electron (muon) has 
a long history. Various relativistic and radiative cor- 
rections have been obtained in analytic or semi-analytic 
form taking into account a deviation of the potential from 
the Coulomb one. Recently interest to an accurate theory 
of the g factor increased. Accurate theoretical approach 
needs to take into account relativistic and radiative cor- 
rections. Here we consider a case when such corrections 
correspond to some modification of a central potential. 
Such a problem is actual because of the Uehling poten- 
tial, which is responsible for a dominant QED correction 
in muonic atoms, finite- nuclear-size effects and some oth- 
ers. The results are presented in a fully relativistic ap- 
proach, i.e. exact in Za for the Coulomb interaction and 
its modifications. 

The methods of calculation of energy levels have been 
successfully developed over the decades. The energy lev- 
els can be calculated with high accuracy by solving cer- 
tain equations analytically or numerically, or using the 
perturbation theory with respect to the corresponding 
correction to the hamiltonian. Calculation of the g fac- 
tor is a more complicated problem. In the traditional ap- 
proach to the calculation of the correction to the g factor 
due to modification of the potential SV, one considers 
both the magnetic field and 5V as a perturbation, thus 
starting with the second order of the perturbation the- 
ory. The numerical approaches are also not easy to apply 
because they usually are much more accurate in determi- 
nation of the energy than the wave function. However, it 
is the latter which is needed to calculate the g factor. 

In this paper, we develop a framework which allows 
one to express the g factor of a Dirac particle bound in 
arbitrary central potential via the binding energy. For a 
small deviation of this potential from the Coulomb one, 
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we derive the general expression for the first correction to 
the g factor due to this deviation. The results are valid 
both for muonic and electronic atoms. 



II. GENERAL RELATIONS 

In the weak homogenous magnetic field B, the cor- 
rection to energy levels of a Dirac particle bound by a 
central potential reads 



AE = ej d 3 r*(r)(7- A)#(r) 
eB ■ (J) 



2m 



where 



x r 



(1) 
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The relativistic units in which K = c = 1 are applied 
throughout the paper; the charge of the bound particle 
is — e, e is a charge of a proton, e 2 = a — 1/137 is the 
fine structure constant. 

Writing the Dirac wave function ip(r) in the form 



A(r)n 

»/ 2 (r)fi 



(3) 



where f2 is the spherical spinor [l| with the angular mo- 
mentum J and orbital momentum L, fl = — (tr • n)fi, we 
obtain for the g factor of a bound Dirac particle 



9 = 



2mn 
+ 1) 



drr z h{r)f 2 {r) 



(4) 



k = j + 1/2 for I > j and k = -(j + 1/2) for / < j. The 
integral in the right-hand side of Eq. (0} can be presented 
as follows (see, e.g., J2|) 



drr 3 fx f 2 = - — 
4m 



1-2/c / drr 2 (f^-fl) 



(5) 
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This formula is valid for the Dirac equation with any 
central potential. Eq. Q has been known for a while 
(see, e.g., Q), but, to the best of our knowledge, it was 
never applied to a non-Coulomb problems. 

Using the identity 



^r 2 (A 2 -/ 2 2 )=( 7 o) 



we have 
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k[1 - 2«( 7 q)] 
2j(j + l) 



(6) 
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Assuming that the potential V(r) is independent of the 
mass of the bound particle, we find 



< 7o ) = /^L 

\ dm 



dE_ 

dm 



(8) 



where H is the Dirac hamiltonian H = 70 (7 -p) +70TO + 
V(r). 

Thus we arrive at an equation, which expresses the g 
factor of the state via its binding energy 



2j(j + 1) 



1 - 2k 



dE 
dm 



(9) 



This equation can be used even in the case when a po- 
tential substantially differs from the Coulomb one, as for 
heavy muonic atoms. In many cases, Eq. @ essentially 
simplifies calculation of the corrections to the g factor of 
a Dirac particle. In particular, it can be applied to such 
problems as finite-nuclear-size effect and vacuum polar- 
ization. 

The equation © becomes essentially simpler if the po- 
tential is close to the Coulomb one, and the deviation 
SV(r) = V{r) — Vc(r) can be treated as a perturbation. 
For the pure Coulomb case, when dEc/dm = Ec/m, 
we immediately obtain the well-known result (see, e.g., 
01) 



Age 



1 - 2k 



7 + n r 
N 



(10) 



2j(j + 1) 
where 

N= V(l + n r ) 2 + (Za) 2 , 7 = V« 2 - {Zaf , 

and n r is the radial quantum number. The correction to 
the energy is 

(11) 



SE= d 6 r^ + (r)8V{r)^{r) 



As follows from the dimensional reasons, the wave func- 
tion in the Coulomb field can be presented in the form 
^(r) = m 3 / 2 $(mr), where ^ is dimensionless. Passing 
to the variable p — mr, we find 



5E= d A p^ + {p)8V{p/m)^{p) 



(12) 



Taking the derivative over m and returning to the vari- 
able r, we obtain 



and 



dSE 
dm 



-— [ d 3 r^ + (r)r 
m J 



dSV(r) 
dr 



tt(r) , (13) 



m + 1) 



2k I d8V(r) 
l + —(r — ^ 
dr 
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Eqs. ©, l|14(l are the basis of our approach. 



III. ANOMALOUS MAGNETIC MOMENT 

It may be interesting to generalize our results to the 
case of a particle with a non-vanishing anomalous mag- 
netic moment. In particular, it is necessary for antipro- 
tonic atoms. The anomalous magnetic moment of the an- 
tiproton is big because of the complicated internal struc- 
ture of the particle. Therefore, it makes sense to consider 
the modification of Eqs. 10, (|14fl due to the anomalous 
magnetic moment separately from the radiative correc- 
tions. 

As known, the modified Dirac hamiltonian for a parti- 
cle with the anomalous magnetic moment in an arbitrary 
electromagnetic field is of the form (see, e.g., 0|) 

H = 10 (j-( P +eA))+ 7o m-eA + 1 -^ joa^F^ , (15) 

where — (1/2)[7 M , 7^], p a is the dimensional anoma- 
lous magnetic moment, F^ v is the tensor of the electro- 
magnetic field; we remind that the charge of the bound 
particle is (— e). 

As follows from Eq. (|15fl . the correction to the energy 
of a bound particle in the homogenous mag netic field © 
is of the form (cf. Eq. Q) 

AE = J d 3 r~${r) {e(~f ■ A) - /i a (£ • B) } <f(r) , (16) 

where T,j = {i/2)ejki(?ki is the spin operator. The wave 
function ^(r) is an eigenfunction of the hamiltonian H15|) 
with — eA (r) — V(r) and Ai = 0. Since the potential 
Ao(r) is spherically symmetric, the wave function still 
has the form ©, and the result for the g factor reads 
(cf. Eq. 0) 



2m 



drr 3 hf 2 



Ma 

2e 



J drr 2 {fl- /|) -2k 
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The radial wave functions /1 and ji satisfy the system 
of equations (cf. [jj]) 



f> + (E-V-m)f 1+ n- J ^-^V^f, = 



f[-(E-V + m)f 2 + 



1 



*L V ' ) h = 



(18) 
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Then we multiply the first equation by r 3 fi and the sec- 
ond one by r 3 fx, sum up the results and take the integral 
over r. Integrating by parts the terms with the deriva- 
tives, we arrive at the following relation (cf. ©) 



drr 3 hf 2 = - — 



1 

4m 



1-2k / drr 2 {f 2 -f 2 ) 



j 'drr 3 V'(f?-f : 



(19) 



Using this formula and Eqs. © and iJSJ, we present the 
result for the g factor in the form 



1 
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+ 1) 



-k(1 + 2a) 



+ (2k 2 + a) 



dE 
dm 



i~/drrV(tf-/ : 



(20) 



where a — — 2m/i a /e is the dimensionless anomalous 
magnetic moment (for a free particle g = 2(1 + a)). This 
result is obtained for an arbitrary potential Aq and is 
exact in the parameter a. We emphasize that the radial 
wave functions f\ and fi and the binding energy E de- 
pend on the anomalous magnetic moment because of the 
equations (|18fl . 

In the nonrelativistic approximation, Eq. I|20|) turns 
into 



1 



(1 - 2k) (a - k) 



" 2j(j + 1) 

- (2 K 2 + 2 K a + a)-!- / drr 3 V f 2 r \ , (21) 
2m 7 J 

where / nr is the radial part of the nonrelativistic wave 
function. This formula is valid even if a ~ 1. 

It is interesting also to consider an expansion of the g 
factor H20JI in the parameter a for arbitrary field strength 
(when the nonrelativistic approximation is not valid). 
The linear in this parameter term reads 



5g a 



m + 1) 

2k 2 d 
m dm 



-2k- 



dE 
dm 



J drr 2 V hf 2 



K 

m 



drr 3 V (f 2 -f< 



(22) 



Here we used the formula for the linear in \i a correction 
to the energy, 



SE n 



2p, a 



drr 2 V hf 2 



(23) 



and took into account that dfj, a /dm = 0. Let us con- 
sider Ea.l)22[l for the pure Coulomb potential. Strictly 
speaking, the equations (|18fl have no sense for a pure 



Coulomb field because of the terms cx 1/r 2 , which lead 
to the phenomenon of falling to the center. It is a conse- 
quence of the point-like source of the field and is absent 
if finite nuclear size is taken into account. However, for 
large quantum numbers which are mostly interesting for 
the experiments, the correction due to the finite nuclear 
size does not change essentially the result obtained from 
Ea. Q22[l for a pure Coulomb field. Using the radial matrix 
element from Ref.H, we obtain 



5g a = 



-2k- 



7 + n r n{Zay 



2j(j + l) I N 
4K 2 (Za) 4 [2K{j + n r ) - N] 
7 (4 7 2 - l)iV 4 



TV 2 



(24) 



where N and 7 are defined after Eq. (|1(J|) . We see that 
the correction to the g factor due to the anomalous mag- 
netic moment has essentially more complicated depen- 
dence on quantum numbers than the leading term I|1U|) . 



IV. FINITE-NUCLEAR-SIZE EFFECT 

As an illustration of efficiency of our method, let us 
consider a finite-nuclear-size correction to the bound- 
electron g factor. For Is state, this problem was solved 
analitically in Ref.[l| in the nonrelativistic approximation 
(Za -C 1), and numerically in Ref.[(| for arbitrary value 
of the parameter Za. In Ref.Q, the result was obtained 
analytically in the next-to-leading approximation in the 
parameter Za, and for arbitrary states. The result of 
Ref.0 allows one to describe well the correction to the 
g factor up to Z — 20. On the other hand, in Ref. 8] 
the approximate formulas for the finite-nuclear-size cor- 
rection to the energy levels, SEf ns , were obtained, which 
gives the result with a relative error of about 0.2% up to 
Z = 100. Using Eq.© and known dependence of Ef ns 
on m, see H, we obtain 



dE fns _ k 2 (2 7 +1) E f , 



3 (j + !) 9m 



j(J + 1) m 



(25) 



Moreover, using our method we can calculate the radia- 
tive correction to Sgf ns for S1/2 state andpi/2 state (and 
arbitrary n r ) coming from the effect of vacuum polar- 
ization. For these states, the finite-nuclear-size effect is 
the most significant. The radiative correction SEj ns to 
the energy Et na was considered in detail in Ref. • The 
correction due to the vacuum polarization SEj^ can be 
represented in the form 



EfnsA 



VP 



(26) 



where Et ns was obtained in Ref. ||, and the explicit 
analytical form of A VP for S1/2 state and P1/2 state is 
derived in Ref.0. Eq. (gSJ with E fns from Ref.Q and 
A VP from Ref. [|| provides the high accuracy of SEj p s 
up to Z — 100. However, the correction A VP was calcu- 
lated in for the mass of the bound particle m being 



4 



equal to the mass M of the particle in the fermion loop. 
Therefore, in order to use Eq. 10 we should know the 
explicit dependence of A VP on m/M (we can set M = m 
only after the differentiation over m) . The corresponding 
expression for A VP can be easily obtained following the 
derivation of A VP in For Za ~ 1, the main contri- 
bution to A vp is given by the logarithmically enhanced 
term 

A VP^MZal ln 2 {mRh (2?) 



where 71 = y/l — (Za) 2 , and R is the nuclear radius. 
This term is the same for S1/2 and P1/2 states and comes 
from the distances l/m, 1/M » r » Ji. Therefore, it 
contains only the logarithmic dependence on the mass of 
the bound particle, in contrast to the power-like depen- 
dence of Ef ns . Within the logarithmic accuracy, we have 



^6EjZ~A VP £.E fns = (2 11 + l)-^. (28) 
dm J am m 

Thus, for Za ~ 1, we obtain 

A b~F VP 

8 9 ^ s = \{2 ll + l)-^. (29) 
1 6 m 

For Za -C 1 the dependence of A VP on m/M is more 
complicated and will be considered elsewhere. 

V. CONCLUSION 

We present a new effective approach to obtain various 
corrections to the bound particle g factor using the cor- 
responding corrections to the energy levels. It should be 
noted that this method is appropriate for the perturba- 
tions of potential-like type. There are also corrections 



which cannot be derived via this approach. The lead- 
ing part of the non-potential contribution is due to the 
virtual Delbruck scattering and it was calculated in Ref. 

As other applications of our method, we can suggest 
the Uehling potential contribution to the g factor of 
bound electrons and muons. In muonic atoms it is one 
of the dominant effects, while in electronic hydrogen-like 
atoms the accuracy is essentially higher (see, e.g., |ll|). 

For muonic atoms, the Uehling-potential correction to 
the energy levels was calculated in the nonrelativistic ap- 
proximation for certain levels in [l2T | . The result exact in 
Za was obtained in for ground state of hydrogen-like 
atom and agrees with Eq. 0. The methods developed 
there can be easily applied to arbitrary levels and thus, 
using Eq. © , one can find the correction to the g factor 
due to the Uehling potential. 

Using the known corrections to the energy levels com- 
ing from the high-order terms of vacuum polarization, we 
can obtain the corresponding contribution to the g factor. 
One correction comes from the Wichmann-Kroll poten- 
tial (which accounts for the higher-order in Za terms in 
the induced charge density) , see Review [lj] . Others are 
the second-order correction with respect to the Uehling 
potential, and first-order correction with respect to the 
radiative correction to the Uehling potential (of order of 
a). 
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